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Abstract. Compressible vortex sheets are fundamental waves, along with shock and rarefaction waves, in 
entropy solutions to the multidimensional hyperbolic systems of conservation laws; and understanding the 
behavior of compressible vortex sheets is an important step towards our full understanding of fluid motions 
and the behavior of entropy solutions. For the Euler equations in two-dimensional gas dynamics, the classical 
linearized stability analysis on compressible vortex sheets predicts stability when the Mach number M > y/2 
and instability when M < \/2; and Artola-Majda's analysis reveals that the nonlinear instability may occur 
if planar vortex sheets are perturbed by highly oscillatory waves even when M > \/2. For the Euler equations 
in three-dimensions, every compressible vortex sheet is violently unstable and this violent instability is the 
analogue of the Kclvin-Hclmholtz instability for incompressible fluids. The purpose of this paper is to 
understand whether compressible vortex sheets in three dimensions, which are unstable in the regime of 
pure gas dynamics, become stable under the magnetic effect in three-dimensional magnetohydrodynamics 
(MHD). One of the main features is that the stability problem is equivalent to a free boundary problem 
whose free boundary is a characteristic surface, which is more delicate than noncharacteristic free boundary 
problems. Another feature is that the linearized problem for current- vortex sheets in MHD does not meet the 
uniform Kreiss-Lopatinskii condition. These features cause additional analytical difficulties and especially 
prevent a direct use of the standard Picard iteration to the nonlinear problem. In this paper, we develop 
a nonlinear approach to deal with these difficulties in three-dimensional MHD. We first carefully formulate 
the linearized problem for the current-vortex sheets to show rigorously that the magnetic effect makes the 
problem weakly stable and establish energy estimates, especially high-order energy estimates, in terms of the 
nonhomogeneous terms and variable coefficients without loss of the order. Then we exploit these results to 
develop a suitable iteration scheme of Nash-Moser-Hormander type and establish its convergence, which leads 
to the existence and stability of compressible current- vortex sheets, locally in time, in the three-dimensional 
MHD. 



1. Introduction 

We are concerned with the existence and stability of compressible current-vortex sheets in three-dimensional 
magnetohydrodynamics (MHD). The motion of inviscid MHD fluids is governed by the following system: 



a t p + v ■ ( P v) = o, 

V-H& m + viw 1\h\ 2 ) = n. 

(i.i) 



dt(pv) + V • (pv ® v - H <g> H) + V(p + \\H\ 2 ) = 0, 



d t H - V x (v x H) = 0, 
ft(p(|(M 2 + lJ f) + e)) + V ■ {pv(±\vf +e+Z)+Hx(vxH))=0 



and 



V-if = 0, (1.2) 

where p, v = (vi, v%, V3), H = (Hi, H2, H3), and p are the density, velocity, magnetic field, and pressure, 
respectively; e = e(p, S) is the internal energy; and S is the entropy. The Gibbs relation 

9dS = de — ^dp 
P 

implies the constitutive relations: 

p = p 2 e p (p,S), 9 = es(p,S), 
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where 9 is the temperature. The quantity c = y/p p {p, S) is the sonic speed of the fluid. 
For smooth solutions, the equations in are equivalent to 



{{d t + V ■ V)p + pc 2 V ■ v = 0, 
p(d t + v ■ V)v + X7p - (V x H ) x H = 0, 
(dt+v- X7)H — (H • V)« + HV -v = 0, 
(d t +v- V)5 = 0, 



(1.3) 



which can be written as a 8 x 8 symmetric hyperbolic system for U = (jp,v,H, S) of the form: 



3 



B (U)d t U + Bj{U)d Xj U = 0. 



(1.4) 



3=1 



Compressible vortex sheets occur ubiquitously in nature (cf. [3*1 l4l ll(JU13l 1211 123U24| ) and are fundamental 
waves in entropy solutions to the multidimensional hyperbolic systems of conservation laws. For example, the 
vortex sheets are one of the core waves in the Mach reflection configurations when a planar shock hits a wedge 
and, more generally, in the two-dimensional Riemann solutions, which are building blocks of general entropy 
solutions to the Euler equations in gas dynamics. Therefore, understanding the existence and stability of 
compressible vortex sheets is an important step towards our full understanding of fluid motions and the 
behavior of entropy solutions to the multidimensional hyperbolic systems of conservation laws, along with 
the existence and stability of shock and rarefaction waves (cf. Majda [T|"], Glimm-Majda [T"""|, Alinhac [2], 
and Metivier [T""j; also see [TTllrS] '). 

For the Euler equations in two-dimensional gas dynamics, the classical linearized stability analysis on 
compressible vortex sheets predicts stability when the Mach number M > y/2 since there are no growing 
modes in this case, and instability when M < y/2 since there are exponentially exploding modes of instability 
(see Miles [""II] for the definitive treatment). The local nonlinear stability for the two-dimensional case with 
M > y/2 was recently established by Coulombel-Secchi [HI E] when the initial data is a small perturbation 
of a planar vortex sheet. In a series of papers, Artola and Majda 3. studied the stability of vortex sheets 
by using the argument of nonlinear geometric optics to analyze the interaction between vortex sheets and 
highly oscillatory waves. They first observed the generation of three distinct families of kink modes traveling 
along the slip stream bracketed by shocks and rarefaction waves and provided a detailed explanation of the 
instability of supersonic vortex sheets even when M > y2- Therefore, one can not expect the nonlinear 
stability globally for two-dimensional compressible vortex sheets. 

For the Euler equations in three-dimensional gas dynamics, the situation is even more complicated. In 
fact, it is well known that every compressible vortex sheet is violently unstable and this violent instability is 
the analogue of the Kelvin-Helmholtz instability for incompressible fluids. 

The purpose of this paper is to understand whether compressible vortex sheets in three dimensions, which 
are unstable in the regime of pure gas dynamics, become stable under the magnetic effect, that is, the 
nonlinear stability of current- vortex sheets in three-dimensional MHD. One of the main features is that the 
stability problem is equivalent to a free boundary problem whose free boundary is a characteristic surface, 
which is more delicate than the noncharacteristic free boundary problems that are endowed with a strict 
jump of the normal velocity on the free boundary (cf. |18l 1111 IB?] : also ["3 EI])- Another feature is that the 
linearized problem for current- vortex sheets in MHD does not meet the uniform Kreiss-Lopatinskii condition 
|17l I18| : see the analysis in Blokhin-Trakhinin [""] and Trakhinin 22 . These features cause additional 
analytical difficulties and, in particular, they prevent a direct use of the standard Picard iteration to prove 
the existence of solutions to the nonlinear problem. 

In this paper, we develop a nonlinear approach to deal with these difficulties in three-dimensional MHD, 
based on the previous results mentioned above. We first carefully formulate the linearized problem for 
the current- vortex sheet to show rigorously that, although any compressible vortex sheet may be violently 
unstable in the regime of pure gas dynamics and does not meet the uniform Kreiss-Lopatinskii condition, 
the magnetic effect makes the linearized problem weakly stable as observed in [51 122|. In particular, we 
successfully establish high-order energy estimates of the solutions for the linearized problem in terms of the 
nonhomogeneous terms and variable coefficients without loss of the order. Then we exploit these results to 
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develop a suitable iteration scheme of Nash-Moser-Hormander type and establish its convergence, which leads 
to the existence and stability of compressible current-vortex sheets, locally in time, in the three-dimensional 
MHD. 

We remark that, in order to establish the energy estimates, especially high-order energy estimates, of 
solutions to the linearized problem derived from the current-vortex sheet, we have identified a well-structured 
decoupled formulation so that the linear problem (|3.1|) - (|3.3|) is decoupled into one standard initial-boundary 
value problem i|3.10[l for a symmetric hyperbolic equations and another problem l|3.11[l for an ordinary 
differential equation for the front. This decoupled formulation is much more convenient and simpler than 
that in (58)-(60) given in [23 an d is essential for us to establish the desired high-order energy estimates of 
solutions, which is one of the key ingredients for developing the suitable iteration scheme of Nash-Moser- 
Hormander type that converges. 

This paper is organized as follows. In Section 2, we first set up the current- vortex sheet problem as a free 
boundary problem, and then we reformulate this problem into a fixed initial-boundary value problem and 
state the main theorem of this paper. In Sections 3-4, we first carefully formulate the linearized problem 
for compressible current-vortex sheets and identify the decoupled formulation; and then we make energy 
estimates, especially high-order energy estimates, of the solutions in terms of the nonhomogeneous terms 
and variable coefficients without loss of the order. Then, in Section 5, we analyze the compatibility conditions 
and construct the zero-th order approximate solutions for the nonlinear problem, which is a basis for our 
iteration scheme. In Section 6, we develop a suitable iteration scheme of Nash-Moser-Hormander type. In 
Section 7, we establish the convergence of the iteration scheme towards a compressible current-vortex sheet. 
Finally, in Section 8, we complete several necessary estimates of the iteration scheme used in Section 7. 

2. Current- Vortex Sheets and Main Theorem 

In this section, we first set up the current-vortex sheet problem as a free boundary problem and then 
reformulate this problem into a fixed initial-boundary value problem, and finally we state the main theorem 
of this paper. 

Let a piecewise smooth function U(t, x) be an entropy solution to Ijl.lll with the form: 

y, t v = f U + (t,x) when Xi > ip(t,x 2 ,x 3 ), 
)U~(t, x) when x\ < ip(t, x 2l x 3 ). 

Then, on the discontinuity front T = {xi = ip(t, x 2 , £3)}, U(t,x) must satisfy the Rankine-Hugoniot condi- 
tions: 

' [m N ] = 0, [H N ] = 0, 
m N [v N ] + [q}= 0, m N [v T ] = H N [H T ], 

m N [^]=H N [v k ] for k= 1,2,3, ( ' } 

m N [e + i(M 2 + ^f)] + [qv N -H N {H- «)] = 0, 

of the function a on the front T, (vn,v t ) (resp. ( 
resp. H) on T, i.e., 

v N := vi - ip X2 v 2 - ipx 3 v 3 , v T = {v Tl ,v T2 ) T := (ip X2 Vi + v 2 ,ip X3 vi + v 3 ) T , 

H N := H x - yj X2 H 2 - yj X3 H 3 , H T = (H Tl , H T2 ) T := {ip x2 H x + H^ Xz II x + H 3 ) T , 

\H\ 2 

rriN — p(vjy — ipt) is the mass transfer flux, and q = p + is the total pressure. 

If rriN = on r, then (C/ ± ,r) is a contact discontinuity for \\.\\ . In this paper, we focus on the case: 

H+ = H N = 0, (2.3) 
H+ WH-, (2.4) 

for which U(t, x) is called a compressible current-vortex sheet. Then the current-vortex sheet is determined 
by (E2Jl and 

& = b + ^r-] = o, (2.5) 



where [a] denotes the jump of the function a on the front T, (vn, v t ) (resp. (-Hat, H t )) are the normal and 
tangential components of v (resp. H) on T, i.e., 



on which generically 

([p],[v T ],[S})^0. (2.6) 

When H = 0, U (t, x) becomes a classical compressible vortex sheet in fluid mechanics without magnetic 
effect. 

Since = on T = {x\ = -0(i, X2, £3)}, it is easy to see that condition (|2.4|l is equivalent to the 
condition: 



H 2 \w( H 2 



on T. (2.7) 



Indeed, when = H N on T, we have 



= ((1 + 4> A X2 )Hf + iP X2 ip x3 Hf^ X2 ^ X3 Hf + (1 + rfejHf) 

which implies that aH+ + (3H~ = if and only if aH^ + [3H^ = hold for k = 2, 3. 

Furthermore, system has the following property: Let (U^ , tp) be the current-vortex sheet to <|1.1|) for 
t e [0,T) for some T > 0. Then, ifV ■ #±(0,35) = 0, we have 

V-H ± (t,x) = for allt e [0,T). (2.8) 

This can be seen as follows: From we find that, on both sides of T, 

(d t + v ± ■ V)(V • iJ ± ) + (V • w ± )(V • iJ ± ) = 0. (2.9) 

On r, we use ipt — v n t° obtain 

d t + v ± ■ V = d t + (v^ + ip X2 Vz + i/} X3 vf)d Xl + vfd X2 + vfd X3 = r ± • (dt, V), 

where the vectors = (l,f^ + ipx 2 v 2 + ^x-^vf ,vf) T are orthogonal to the space-time normal vector 
n = (ipt> — 1, ^Pxz) T on r. Thus, ^ + w ± • V are tangential derivative operators to T. Thus, equations in 
(12.91) are transport equations for V • H on both sides of Y, which implies l|2.8|l . 



Set 

D(X, U) := ( D{X Q U) J ) with D(X, U) := I A P ff ^J 3 -pAJ 3 

\ 0,3x1 —A/3 I3 

As in [22, property (|2.8|l implies that system (|l.l(l - (|1.2(l is equivalent to the following system on both sides 
ofT: 

3 

D(X ± ,U ± )(B (U ± )d t U ± + Y,B J (U ± )d X] U ± ^j +\ ± G ± V-H ± = 0, (2.10) 

3=1 

provided V • ff ± (0, x) = 0, where 

G ± = -(l,0,0,0,i/ ± ,0) T , 
and A ± = A ± (w 2 t , vf, H^, H^) will be determined later so that 

(^) 2 < „ + ,_if?' +l2 (2- 11 ) 



are satisfied for the sound speeds c^. 
We rewrite system (12. 10f) as 



3 

A^U^dt^ +Y J MU ± )d Xj U ± = 0, (2.12) 
j'=i 

which is also symmetric hyperbolic. Then the problem of existence and stability of current-vortex sheets of 
the form 1)2.1(1 can be formulated as the following free boundary problem: 
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U\ t =o = 



Free Boundary Problem: Determine U^(t,x) and a free boundary T — {xi — ^(t,X2,x 3 )} for t > 
such that 

'MU ± W ± + Y? J=1 MU ± )9* i U ± =0 for ±( Xl -iP{t,x 2 ,x 3 ))>0, 

Uq(x) when x\ > ipo{x2, x 3 ), (2-13) 
Uq(x) when x 1 < ipo(x 2 ,x 3 ) 

satisfying the jump conditions on T: 

^t = 4, H± = 0, [p+^L]=0, ff+^ff-, (2.14) 

where i/jq = tj)(0, X\, x 2 ), and condition (I2.11|) is satisfied for \^ given by <|3.9fl later. 

For a piecewise smooth solution U(t,x) of problem (|2 . 1 Hf) — (|2 . 14p with smooth (U^ji/j) for t £ [0, T), it 
can be similarly verified that, if V • i? ± (0, x) = 0, then V • H (t, x) = for all £ € [0, T). Hence the equations 
in Q2.13|l and (|1.1JI - H1.2|I are equivalent for such a solution. 

Furthermore, unlike the shock case, the new difficulty here is that the free boundary T in problem (|2.13|) - 
(I2.14() is now a characteristic surface of system (|2.13|) . rather than the noncharacteristic boundary surface 
endowed with a strict jump of the normal velocity as in the shock case (see P ITU1 ITT1 IT%1 H*^] V 

To deal with such a free boundary problem, it is convenient to employ the standard partial hodograph 
transformation: 

U = t, X 2 =X 2 , X 3 =X 3 , ^ 
[Xi = ^ ± (t,Xl,X2,X 3 ) 

with , J ,± satisfying 

i±d xi ^ ± > k > 0, ,„ „, 

I - 2.16) 

| k *+|s 1 =o = *~k=o = ^(t,x 2 ,x 3 ) 

for some constant k > 0. Under H2.15[l . the domains f2 ± = {±(xi — ip(t, x 2 , x 3 )) > 0} are transformed into 
{xi > 0} and the free boundary T into the fixed boundary {x\ — 0}. 

Then we define C/ ± (i, x) := J7 ± (t,x). From the first jump condition in (|2.14|) . the natural candidates of 
\l/ ± are those that satisfy the following eikonal equations: 

5 t * ± - vf +v 2 t d X2 ^ ± +vfd X3 V ± = in{xi>0}, (2.17) 

where we have dropped the tildes in the formula for simplicity of notations. 
It is easy to check that, under l|2.15|l . U ± (t,x) = U ± {t 1 x) satisfy 

L([/ ± ,* ± )J7 ± = in {£!>()}, (2.18) 

the boundary conditions on {x\ = 0}: 

f* + |xi=0 = *"|x 1= = V 1 , 

\B(U+,U-, ^)| B1=0 = 0, 
and the initial condition at {t = 0}: 

(tf±V0lt=o = (^(x), Vo(x 2 ,x 3 )), (2.20) 
where the tildes have also been dropped for simplicity of notations, 

3 

L(U, *)V = A (U)d t V + A X (U, *)d Xl V + A,j{U)d X] V 

with Ax(U, *) = ^f(Ai(U) - d t ^A (U) - E - =2 d x ^MU)), 

B(U+, U-,1>) = (dtil> - U± N , U± N , q+ - q~y 
with U± N = TJf - ^ X2 Uf - ip X3 Uf, U± N = Uf - i) X2 Uf - ip X3 U±, and 



(2.19) 



q = Ut+ ^\U H \ 2 for U H = {U 5 , U e , U 7 ) 
5 



a(xi) 



HU, T :=( J2 W Ma 9 k Xl n\\h ( a T) y 



With these, the free boundary problem has been reduced into the fixed initial-boundary value problem 

To solve ipT7)l - ip3D|l . as in PJ0|n|, it is natural to introduce the tangential vector M = (Mi, M 2 , M 3 ) 
of {2:1 = 0}: 

Mi = o-fci)^! M 2 = d X2 , M 3 =d X3 , 

with 

a;i when < xi < 1, 
2 when X\ > 2, 

and the weighted Sobolev spaces defined on Qj- := {(t, x) e [0, T] x 1R 3 : x\ > 0}: 

S*(nr) := {« G X 2 (^t) : e~ fJ ' t M a d^u £ L 2 (n T ) for |a| + 2fc < s} 
for all s G W and > 0, imposed with the norms 

nut :=( A^oim 1 ' 2 . 

Jo 

where 

IKV)II2, M := E M 3( - |o| - 2 * ) l|e-' ,t M '^ 1 tt(t,.)||i a . 

|a|+2/c<s 

We will also use similar notation as above for the spaces with /i = 0, B s (flT), whose norm is defined by 

1 1/2 

1^ u x 1 u \\L^(n T ), 

a|+2fc<s 

Also denote 6ft T := {(t,x 2 ,x 3 ) : t G [0, T], [x 2 ,x 3 ) G J? 2 }. 

Then, when the initial data (U ,ipo) is a small perturbation of a plane current-vortex sheet (i7 ± ,'0) for 
constant states and -0 = satisfying (|2.3|) - l|2.5|) . we have the following main theorem of this paper. 

Theorem 2.1 (Main Theorem). Assume that, for any fixed a > 14 and s £ [a + 5, 2a — 9], ifte initial 
data functions ipo G H 2s+3 (IR 2 ) and Uq — G _B 2 ^ s+2 ^(5?^_) satisfy the compatibility conditions of problem 
(|2.17() - (I2.20() up to order s + 2. 77ien iftere exists a solution (U^, \t' ± ) 0/ i/ie initial-boundary value problem 
$FT[\ - RFM such that 

(C/± - tf± - ¥±) G B Q (ft T ) arid ^ = $±1^=0 G H a -\bn T ) 

for some functions and satisfying 

U± - U ± , ^Tn€ B S+3 (]R + x JR^. ). 

We remark that, since ±<9j 1 >I' ± > k > 0, the corresponding vector function of the solution (L/ ± , , & ± ) 
in Theorem 2.1 under the inverse of the partial hodograph transform is a solution of the free boundary 
problem, i.e. the current-vortex sheet problem, which implies the existence and stability of compressible 
current-vortex sheets to l|2.1|l - (|2.2|) under the initial perturbation. To achieve this, we start with the linear 
stability in Sections 3-4 and then establish the nonlinear stability by developing nonlinear techniques in 
Sections 5-7. Some estimates used in Sections 5-7 for the iteration scheme of Nash-Moser-Hormandcr type 
are given in Section 8. 

3. Linear Stability I: Linearized Problem 

To study the linear stability of current-vortex sheets, we first derive a linearized problem from the nonlinear 
problem H2.17fl - H2.20fl . By a direct calculation, we have 

4- (L(U + sV, * + a$)(ff + sV)) \ s=0 = L(U, ¥)W + E(U, *)W + (L(U, tf)Z7), 

where 

W = V ~^ U 
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is the good unknown as introduced in |2] (see also ^2 US]) an d 



Then we obtain the following linearized problem of (|2.17|) - l|2.20|) : 

L{U ± ,^ ± )W ± + E{U ± ,^ ± )W ± = F ± in {xi > 0} 
with the boundary conditions on {x\ = 0}: 

'4> t ~ {Wt - i> X2 w+ - fl> Xa Wt) + vtt*, + ufK, = hf, 
wt - 1>**wt - 4>* a w? - Uf<t> X2 - u±<t> X3 = hf, 
w? - wr + T,U(u^w+ ~ ujwj) = h 3 , 

and the initial condition at {t — 0}: 

(^ ± »| t=0 = 0, 
for some functions and h := (hf, hf,h^) T . 

To separate the characteristic and noncharacteristic components of the unknown W , we introduce = 
J(C/ ± , \E ,± ) as a 8 x 8 regular matrix such that 

X ± = {J ± )- 1 W ± (3.4) 

satisfies 

'xf^wf + Y^Uufwf, 
xf = wf - y±wi ~ vf 3 w£, 



(3.1) 



(3.2) 



(3.3) 



^ ■2 



XT = WT - <f ± Wt - * ± Wt 



. , x± , xf , xf ,xf) = (wt , wt , wt , wt ,wt). 

Then, under transformation (|3.4II . problem H3.1fl - I|3.3|l for (W^ , <j>) is equivalent to the following problem for 
(X ± A): 

X(J7 ± , * ± )X ± + E(U ± ,^ ± )X ± = F ± in {xi > 0} 
with the boundary conditions on {x% = 0}: 



(3.5) 



b t - Xf + Uf^ + Ut<f> X3 = hf, 



u 6 4> X2 - u 7 4> X3 = h 



± 

2 • 



Xf ~ X 1 — /13, 



and the initial condition at {t = 0}: 
where F ± = (J ± ) T F ± , 



(X ± ,<j>)\ t=o = 0, 



(3.6) 



(3.7) 



i(t/ ± , *±) = Aoiu*, *±)a t + * ± )^ J 

i=i 

with A l (U ± ,^ ± ) = (J ± ) T A 1 (C/ ± ,* ± )J ± , A^l/*,**) = (J ± ) T A 3 (J7 ± )J ± , j ^ 1, and 

£7(17* * ± )X ± = ( J ± ) T J B([/ ± , * ± ) J ± X ± + (J ± ) T (L(U ± ,V ± )J ± )X ± . 

For simplicity of notations, we will drop the tildes in (|3.5(l . By a direct calculation, we see that Ai(J7 ± ,* ± ) 
can be decomposed into three parts: 

A 1 (J7 ± ,* ± ) = Af 

with 



r I A :.,0 , ^±,1 , ^±,2 



±,0 



1 



a 

a T 7x7 
r± 



.4 



±,2 
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where a = (1, 0, 0, -A±, 0, 0, 0), U* N = Uf - V X2 Uf - ^ x ,Uf, and Ufj N - „ B 



d Xl V± 1 



When the state (£/ ± ,\£ ,± ) satisfies the boundary conditions given in i|2.19[l . we know from (|3.5[) that, on 



Ai(U-,$-) J \ X~ J ' \ X 

when [Xi] = on {x\ = 0}. 

Furthermore, from (|3.6[) . we have 

[X 2 - \X 5 ] = <j> X2 [U 3 - XU 6 ] + <t> X3 [Ui - AIM - [hi + Xh 2 ]. 

Thus, from assumption l|2.4[l that is equivalent to (|2.7[) . there exists a unique A* = A ± (w^ : , vf, H^, H^ 1 ) 
such that 

= A + ($W( 2 2 - V (3-9) 



\ H 3 J V #3 / 

that is, [1^3 — Al^] = [C4 — AC/7] = 0, which implies 

[X 2 -XX 5 ] = -[h 1 +Xh 2 ]- 

Therefore, with the choice of X ± in 13. 9|) . problem (|3.5|) (|3.7|) is decomposed into 

' L(U ± ^ ± )X ± + E(U ± ,^ ± )X ± = F ± in {x x > 0}, 
[X 2 -AX 5 ] =-[h x + Xh 2 ] on{.T 1= 0}, 
X+ - Xf = /i 3 on {arx = 0}, 



(3.10) 



X ± | t<n = 



which is maximally dissipative in the sense of Lax-Friedrichs ^J^]> and 



" X > - (U >- U ' ] u} u 7 - \x}- + " f °° (X1 = ° h (3-11) 



it=o = 0. 



Remark 3.1. To have the identities in (|3.11(l . it requires the compatibility, that is, the right-hand side 
of 1|3.11|) for the plus sign is equal to the term for the minus sign. This is also guaranteed by the choice of 

a ± in iprg|> . 



4. Linear Stability II: Energy Estimates 

In this section, we establish the energy estimates for the linear problem (|3.1(l - (|3.3|l . which implies the 
existence and uniqueness of solutions. We know from Section 3 that it suffices to study the linear problem 
(I3TTU1) - (I3TTT1) . 

As noted as above, when the state (C/ ± ,4' ± ) satisfies the eikonal equation i|2.17|l in {x\ > 0} and the 
boundary conditions on {x\ = 0} in l|2.19|> . we know that the boundary {x\ = 0} is a uniform characteristic 
of multiplicities 12 of the equations in (|3.10|) . 

First, we have the following elementary properties in the space B s (Q,t), whose proof can be found in 

Lemma 4.1. (i) For any fixed s > 5/2, the identity mapping is a bounded embedding from B s (Qt) to 
L™(n T ); 

(ii) If s > 1 and u 6 B s (Q,t), then u\ Xi —q € iP (6fir)> the classical Sobolev space. Conversely, if 
v G H° (Wlr) for a fixed a > 0, then there exists u £ B a+1 (rtT) such that u\ Xl =a = v. 

Denote by coef(t,x) all the coefficient functions appeared in the equation in l|3.1U|l . and coef(t,x) — 
coef(t,x) — coe/(0). With the decoupled formulation H3.10fl - H3.11f) for the linearized problem, we can now 
employ the Lax-Friedrichs theory ^]^| to establish the desired energy estimates, especially the high-order 
energy estimates, of solutions for the linear problem. 
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Theorem 4.1. For any fixed sq > 17/2, there exist constants Co and liq depending only on \\coef\\ So ,T 
for the coefficient functions in i|3.10[) such that, for any s > sq and /i > Lio, the estimate 

nmxjX ± (t)\\l^ + Li\\X ± \\l^ T 

< jiWFHl^T + \\h\\ 2 H ^ (mT) + ll^/II^Tdl^ll^.T + INI^o +1(6 n T ))) 41 

holds provided that the eikonal equations <|2.17[1 and 

Uh,n ■= Uf ~ d^Uf - d x ^Uf = 

are valid for (U^, \I r± ) in a neighborhood of {x\ = 0}, where h = (hf, h 2 , hs) T and the norms in H^bflx) 
defined as that of B^(£It) with functions independent of x\. 

Proof: The proof is divided into two steps. 

Step 1. We first study problem l|3.10|l with homogeneous boundary conditions: 

hf = hf = h 3 = on {xx = 0}. (4.2) 
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A 2 { A. 



A fx=[l oj' Af^=A^=Af^ = (4.3) 

in a neighborhood of {x\ = 0}. 

Thus, from l|3.1()Jl and l|4.2|l . we find that the vector function 

Y ± = (P ± y 1 X ± 

satisfies the following problem: 

'CQJ^, * ± )F ± + £ ± {U ± , * ± )y ± = F± in {x x > 0}, 

[Y 1 ] = [Y 2 }=0 on{a; 1 =0}, .4.4; 
,y ± |t<o = 0, 

where ^ = (P ± ) T F ± , 

3 

£(t/ ± ,vi/ ± ) = ^ ± a t + ^^ ± ^ j 

3=1 

with Af = (P ± ) T AfP ± ,0< j < 3, and £ ± = (P±) T EP ± + (P±) T £ 3 3 =0 Afd Xj P ± . 

L 2 - estimate: With setup (|4.4() . we can apply the Lax-Friedrichs theory to obtain the following L 2 — estimate: 

m^ll^WII^+^II^II^T^^II^IIw (4.5) 

when fi > fiQ > 0, where 

Co = Co(||coe/||z,~, ||V • coe/||z,oo). 
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Estimates on the tangential derivatives. From 14. 411 . wc find that M a Y ± , \a\ < s, satisfy 

'CM a Y ± + £ ± M a Y ± = M a F± + [C + £ ± ,M a ]Y ± in {x t > 0}, 
[M a yi] = [M a Y 2 ] =0 on {xi = 0}, (4.6) 

K M a Y ± \ t <p = o, 

which, by the Lax-Friedrichs theory again, yields 
ma* r \Y ± (t)\l fl + vlY^l^ 

~ r ( 4 - 7 ) 

< -tG^I^t + (llcoe/l^) 2 !^!^ + HV^III^Icoe/l^ + \d x ^\ 2 8 ^ T ), 

where Yf 1 — (Y^,Y^), the norms | ■ \ Sifl and | ■ \ s ,fi,T are the special cases of the norms || • || S)/J and || ■ \\ s ,^,t 
respectively without the normal derivatives d£ , and 

Hlr= E ll^lU-cnr) + E ll^^ulU-cnr). 

M<2 I"I<1 



Estimates on the normal derivatives. Set Aj — diag[Aj,Aj ],0 < j < 3, with 

_4± _ / A il ^ 



12 

3 I 4±>j 4±J 

^21 ^22 



From |@3Jl-g3J|, we find that Fj = (Y^, Y,f) and *j = (*3 > ■ ■ ■ > Y s f) satisfy the following equations: 

A^{ l d xl Y^ = coef ■ MY ± + coef ■ Y ± + jf, (4.8) 

and 

{V 3 4 ±J r) V ± 4- r ± V ± — T*- — V 3 4 ±J 7) V ± ^r< ± Y ± 
Z^i j=0 ^22 u ^j 1 I ' ^- / 22 1 1 -'J Z^j=0 /1 21 u Xj I I l ~ y 21 I I i (4 9) 

^|*<o = 0, 

where the linear operator X)j=o ^22'^^ + C22 is symmetric hyperbolic and tangential to {xi = 0}. 
Thus, from 1)4. 8|) . we obtain 

\d xl Y±{t)\ s _ lttl < \jf(t)\ s ^ + Wcoefh^T^it)^ + WY^L^lcoefit)^ (4.10) 

with || • |j Li P (T) denoting the usual norm in the space of Lipschitz continuous functions in Ot and, from (|4.9|l . 
we deduce the following estimate on d Xl Y^ by using the classical hyperbolic theory: 

max \d Xl Y±(t)\ 2 s _ 2>fl + »\d Xl Y±\ 2 s _ 2 ^ >T 



< + WcoefWl^dYfU^T + Id^lU,^) (4.11) 

+ {\\coef\\* T \Y±\^ T Y + (ll^HJIcoe/l^) 2 ), 

where 

\ u \s^T = \u\s,n,T + \u xl \ a -2 llit T. 

Combining estimates (|4.7f) with (|4.1(J|) - (|4.11|) . applying an inductive argument on the order of normal 
derivatives 9* , and noting that all the coefficient functions of C and £ ± in problem (|4.4|) depending on 
(U ± , V t , 2; * ± ) and {U ± , V fja; C/ ± , {d^ V"^ 2 a . 3 'S f± } Qfl <i ! | Qf2 |<2, ai +|Q :2 |<2), respectively, we can conclude esti- 
mate (I4.1|) for the case of homogeneous boundary conditions by choosing s > 17/2. 

Step 2. With Step 1, for the case of nonhomogeneous boundary conditions in (|3.1U|) with (hf, h 2 , € 
£P +1 ([0,T] x K 2 ), it suffices to study (l3~TU|) when F ± = 0. 
By Lemma 4.1, there exists Xq g B s+2 (Qt) satisfying 

[X ,2-\X^]=-[h x + \h 2 ], [X 0tl ]=ha on{a;i = 0}. (4.12) 
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Thus, from (13 . 1 L)|) . we know that F ± = X ± — Xq satisfies the following problem: 

'L(U ± ,^ ± )Y ± + E{U ± ,^! ± )Y ± = /± in {xi > 0}, 
[Yx] = [Y 2 - XY 5 ] =0 on {x x = 0}, (4.13) 

with /± = -L(C/ ± , * ± )X ± - E^, * ± )X ± G B S {Q T ). 

Employing estimate (|4. 1|) in the case h = for problem H4.13[l established in Step 1 yields an estimate 
on Y^, which implies the estimates of = Y^ + Xq. Combining this estimate with that for the case of 
homogeneous boundary conditions yields (|4.1|l for the general case. 

By choosing /i> 1 and T so that /iT < 1 in 14.1(1 , we can directly conclude 

Theorem 4.2. For any fixed sq > 17/2, there exists a constant Cq > depending only on ||eoe/|| SOi T 
such that, for any s > sq, the estimate 

11**112,1 < Codl^H^T + WHW( b n T ) + I|c^/II2 i t(II* ,± II2o,t + IWI^ + 1(6 n T ))) (4.14) 

holds. 

5. Nonlinear Stability I: Construction of The Zero-th Order Approximate Solutions 

With the linear estimates in Section 4, we now establish the local existence of current-vortex sheets of 
problem (|2.17I) - (|2.2U|) under the compatibility conditions on the initial data. From now on, we focus on the 
initial data (U^,ipo) that is a small perturbation of a planar current-vortex sheet (U^, ip) for constant states 
U ± and $ = satisfying 

We start with the compatibility conditions on {t — x\ = 0}. For fixed k G IN, it is easy to formulate the 
j*' 1 — order compatibility conditions, < j < k, for the initial data (U^ipo) of problem l|2.17|l - H2.20|) . from 
which we determine the data: 

{&t £/ ± | t =: El =o}o<j<fc and {d 3 t ip\t=o}o< 3 <k+i- (5.1) 
For any fixed integer s > 9/2 and given data ipo £ H S ~ 1 (1R 2 ) and with 

(ft = U± -U ± e B S (R 3 + ), 
by using Lemma 4.1, we can extend ipo to be 'J'J G B s (M 3 i _) satisfying 

which is possible by using assumption l|2.4|l . Let = 4-^ ± x\ be the initial data for solving ^ r± from the 
eikonal equations l|2.17|l . 

Let Uf{x) = dlU^a and $>f = fl^^o for any j > 0. 

For all k < [s/2], from problem l|2.17|) - (|2.2U|) . we determine 

{Uf}i< <k and {^} <j<fc+i ( 5 - 2 ) 

in terms of (Uq , ) and 

[uf e B S -^{R%), 1 < j < k, 

\vf G B S - 1 (M%), G B s - 2 ^-^(]Rl), 2 < j < k + 1. 

Furthermore, we have the estimate 

[a/2] [a/2]+l 

^11^11^ — ^(113)+ H*i = Hs'-2j + 2(iR3 ) + H^flls— 1(1?3) 

< Co(||t / o ± lls s (K+) + II^oHh— 

with C depending only on s, ||{/ ± ||^i,°o(. K 3 ), and ll^ollw 1 ."^ 2 )- 

li 



(5.3) 



(5.4) 



Set si — [s/2] + 1. We now construct the zero-th order approximate solutions (£7^=, vjr^). First, we 
construct 



$± g B S1+1 (]R + x ml) 
satisfying aj?#±| t=0 = «f ,0 < j < s u and *+U 1= o = k=o, and 

U±eB s i(lR. + xlR 3 + ), 1^2, 
satisfying S/f7± | t=0 - E/J„0 < j < si - 1. 



(5.5) 



(5.6) 



Moreover, we have 



' K* ~ d *^tUt# - d *s*t U a,7 = m {x x > 0}, 
[U a ,i + h\Ua, H \ 2 } =0 on{x 1= 0}, 



(5.7) 



where U± = U ± + f/±, ¥± = *± ± x u and E^* = (C/ Q , 5 , U afi , U a , 7 ) 
Finally, we construct U± 2 e B Sl (R + x R 3 + ) by requiring 



{x x > 0}. 



(5.8) 



Therefore, we conclude 



Lemma 5.1. For any fixed integer s > 9/2, assume that the initial data (U^ipo) = {U^ — U ± ,ipo) G 
i3 s (J?+) x £F _1 (jR 2 ) is bounded in the norms. Then there exists (U^^t) such that U± = U± - U^ 1 € 
B[V2]+1(JR + x $r± = ^r± =p Xl e B[ s / 2 ]+ 2 (iR + x JZ* ), 



dt(L(U±,^)U±))\ t=0 = /or < j < [s/2] - 1, 



(5.9) 



and ^oZd. 



With this, we set ip a — Ui=o- From the compatibility conditions of problem (|2.17|) - I|2.2U|) and Lemma 
5.1, we have 



%{B(U+, U-^ a ))\ t=0 = for < j < [s/2] 



(5.10) 



Set 



V ± = U ± - uf, 



(J)* = vp* — 1]/= 



(5.11) 



Then it follows from Lemma 5.1 and 15.10f) that problem (|2.17() - (I2.20() is equivalent to the following fixed 
initial-boundary value problem for (l/ ± ,$ ± ): 



C(y±, ^> ± )V ± = /± in {t > 0, Xl > 0}, 
8(V ± ,$ ± ) = Q in {* > 0, xx > 0}, 

$+U 1=0 = $~L=o = <t>, 

B(V+,V-,<j)) =0 on{a;i = 0}, 
_V r± | t <o = 0, $±| t < =0, 

where /± = -L([/±, *±)C/±, 



(5.12) 



and 



6. Nonlinear Stability II: Iteration Scheme 

From Theorem 4.1, we observe that the high-order energy estimates of solutions for the linearized problem 
(I3.1U|) - (|3.11|) in terms of the nonhomogeneous terms and variable coefficients coef keep the same order. 
Based on these estimates and the structure of the nonlinear system, we now develop a suitable iteration 
scheme of Nash-Moser-Hormander type (cf. f° r our nonlinear problem (]5. 121) . 

To do this, we first recall a standard family of smoothing operators (cf. [21III|) : 

{Sg} e>0 : Bl(n T ) — » n s > B s ^ T ) (6.1) 

satisfying 

'\\Seu\\ s , T < C6^-^+\\u\\^ T for all s,a> 0, 
\S e u-u\\ s>T < C6 a - a \\u\\ a>T for all s £ [0, a], (6.2) 
\-^S e u\\ S;T < C0 a - Q, - 1 ||u|| aiT for all s,a > 0, 
and 

|| {S e u+ - S e u-)\ Xl =o\\ H *(.bn T ) < CO^ 1 -^ \\(u+ - U-)\ Sl= o\\ a , T for all s, a > 0. (6.3) 
Similarly, one has a family of smoothing operators, still denoted by {Sg}e>o acting on H s (bQT) 1 satisfying 
(E3 as well for the norms of H s (bn T ) (cf. ElE]). 

Now we construct the iteration scheme for solving the nonlinear problem Ij5.12|l in M+ x M 3 . 

The Iteration Scheme: Let V ±fi = $ ±,D = 0. Assume that (V ± ' k 1 ^ ± ' k ) have been known for 
k = 0, . . . , n, and satisfy 

{(Y ± ' k ,* ± > k ) = in{t<0}, 
\$+' fc k=o = $-' fc k=o = <t> k on{a:i = 0}. 
Denote the (n + l)*' 1 approximate solutions to l|5.12|l in M + x M 3 by 

y±.n+! = y±." + 6V ± ' n , $±' n+1 = $±» n + 8^'", c/) n+1 = (f> n + 8<p n . (6.5) 



Let 6o > 1 and 9 n = ^JO^+n for any n > 1. Let Sg n be the associated smoothing operators defined as 
above. We now determine the problem of the increments {8V ±,n , <5$ ±,n ) as follows: 

Construction of 8V ziz ' n : First, it is easy to see 

where y ±,n+ ^ are the modified states of Sg n V ± ' n defined in (|6.24|) - (|6.25|l later, which guarantees that the 
boundary {x\ = 0} is the uniform characteristic of constant multiplicity at each iteration step, 

e,(Lr±+y ± ' n +3,*±+s f ,„$±,») ^ ^ 

= L{Ut + V ±,n+ *,*± + S 0n ^ n )8V^ n + E{Uf + V ± - n+ ^,^t + S , e „$ ± '") ( 5I> ± '™ 
are the effective linear operators, 

«5F±.« = *y±.» - s ^ n dx ^ U l +V± ' n+ j\ (6.8) 

are the good unknowns, and 



4 
3 = 1 



e±,n = Y j e% (6.9) 



(6.10) 



are the total errors with the first error resulting from the Newton iteration scheme: 
the second error resulting from the substitution: 
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the third error resulting from the second substitution 

(3) _ jl 

and the remaining error 



«£. = ^s^^s^SV^^") - L> ^ 4 ^ ±1 . ) («V ± -, (6-12) 



e ±- = dx ^r+Se *±,") 9MU " + V± ' nH >*« + S o^ n )(Ut + V^)). (6.13) 
Similarly, we have 

B(V + < n+1 ,V-' n+ \c/> n+1 ) -B(V + ' n ,V- n ,r) = B' (u ± +v± , n+ r ^ +Sg ^ n) (SV + > n ,SV-' n M n ) + S n , (6-14) 
where 

4 

e~n = E e t } ( 6 - 15 ) 

J'=l 

with the first error resulting from the Newton iteration scheme: 

e« =B(U+ + V+' n+1 , U~ + V-' n+ \^ a + - B(U+ + V+' n , U~ + V^ n ^ a + <j> n ) 

- B^y^^iSV^W-^Sr), (6 - 16) 
the second error resulting from the substitution: 

el 2) = B[^ +v ^ >Mn) (SV + ' n ,SV-' n ,8r) ~ B^ +s ^ v± .^ m+8 ^ ) {W + ' n ,6V-»S4r) i (6.17) 
the third errors resulting from the second substitution: 

and the remaining error: 

Using 1)6. 6fl and 16.14|l and noting that 

£(F ±, 0j § ±,oj V ±,o = 0) B (y+> , v-'°, t ) = 0, 

we obtain 



(6.20) 



Observe that, if the limit of (U ±,n , $ ±,n ) exists which is expected to be a solution to problem l|5.12[l . then 
the left-hand side in the first equation of (|6.2UI) should converge to , and the left-hand side in the second 
one of H6.20f> goes to zero when n — + oo. Thus, from H6.20|) . it suffices to study the following problem: 

_SSV +,n >SV~' n M n )=9n onM2 T , (6.21) 



((7^+y ± '" + 2,v a +s () „0™) 
V r± ' n |t<o = 0, 

where and # n are denned by 

n 7i — l n n—1 

E/f +^(E e ±.^ = ^„/ a ± , E*+^(E e ^) = > ( 6 - 22 ) 

J=0 j=0 j=0 j=0 

by induction on n, with /q = Sg f^ and 50 = 0. 

We now define the modified state V ±,n+ ^ to guarantee that the boundary {x\ = 0} is the uniform 
characteristic of constant multiplicity at each iteration step ()6.21(l . To achieve this, we require 

(B(V+' n+ K Sg„ <t> n ))f \ Xl =o = fori = 1,2, and all neW, (6.23) 

14 



which leads to define 
and 



Vp n+h = Se n V±' n fori e {3,4,6,7}, 



v ±,n+i = ^g^i.n) + q x ^± + Sf} ^±^ v ±,n+1 + g ^± + S9w$ ±,„) K ±.«+J 
V 5 ±,n+i = fl xa (*± + Sfl B $ ± ' n )V r 6 :t ' B+ * + d X3 (*t + S e „$±< n )V T ±>n+i 



(6.24) 



(6.25) 



Construction o/ J^^'" satisfying <5<I? '"| Xl= o = 8(f) n . Clearly, we have 



£(V ± > n +\$ ± < n+1 )-£(V ± < n ,$ ± > n )=£' .! (5V r± '",,5$ ± '") + e ± 



where 



,3 ^ r 3 



is the linearized operator of £ and 



with the first error resulting from the Newton iteration scheme: 



, ±i „ = £ (F ± ' n+1 , $±>" +1 ) - £ (F ±,n , $±' n ) - £ ( V±.»,*±.»)(*V ±,n , <5$ d 



(6.26) 
(6.27) 

(6.28) 
(6.29) 



the second error resulting from the substitution: 

4 2 ,l=^(V±.^*±.»)(^ ± '^^ ± ' n )-^„y±.^ Se „*±.»)(^ ± ' n .^ ± • 7^ ) 
=^ 2 ((7 - &J* ±,n )^3 ±,B + d X3 ((I - Se n )^' n )SVt' n 

+ V~ S r fl„)V 3 ±,n «x a («* ±,n ) + (/ - ^„)V 4 ± ' n fl x ,(J$ ± '") J 

the third error resulting from the second substitution: 

which vanishes due to (|6.24|) and that all the coefficients of the linearized operator £[ v ± are independent 
of (V 2 1 , V5 ), and the remaining error: 



(6.30) 



(6.31) 



,(4) _<•/ 



{5V^ n ,5^ n )-£' 



(jy ± '",(5$ ± ^). 



Thus, we have 



which leads to define (5$ ±,ra that is governed by the following problem: 



v( J$ ± ' n | t <o = 0, 



where are defined by 



n-l 



(6.32) 



(6.33) 



(6.34) 



(6.35) 



fe=0 



fc=0 



by induction on n starting with h Q = 0. 
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By comparing (|6.34|) with H6.21[l . it is easy to verify 



(6.36) 



From H5.9|) - (|5.1U|I and (|6.34|) . we know that the compatibility conditions hold for the initial-boundary 
value problem (|6.21[1 for all n > 0. 

The steps for determining (SV ± ' n , 5$ ±! ") are to solve first 5V ± > n from (|6.21|) and then (5$ ± '™ from Q6.34p . 
Finally, SV ,n can be obtained from l|6.8[) . 

7. Nonlinear Stability III: Convergence of the Iteration Scheme 
and Existence of the Current- Vortex Sheet 

Fix any so > 9, a > sq + 5, and Si G [a + 5, 2a — So]. Let the zero-th order approximate solutions for the 
initial data (U^ipo) constructed in §5 satisfy 

||t^||. 1+ 3,T + ||*J||. 1+ 3,T + ||/^|| il -4,T<* ) \\f±\\ a+h T/6 is small (7.1) 

for some small constant S > 0. 

Before we prove the convergence of the iteration scheme <|6.21|) and (|6.34|) , we first introduce the following 
lemmas. 



Lemma 7.1. For any k £ [0,n], 

(T/ ±,fc 5 $±,fe)|| s T + ||0 fe || ffs - 1(bf2T) < (7(56* 



(s-q) 4 



for s S [s , si], s 7^ a, 



^llff°-i(6J2 T ) < C5\og6 k , 



< C56» 



(s-a)-| 



for s > sq,s ^ a, 



(7.2) 



| Sg h *±> k )\\ B ,T ■ 

|(5 efc ^ ± ^',^ fc $ ± ' fc )|U,T + ||S flfc ^|| ff «-i (6nr ) < CTlog^, 

I ((J — S'e l< .)F ±,fe , (J — Sg k )^ ± - k )\\ Sy T + — Sg k )(f> k \\ H s-i( mT ) < C56 s k ~ a for«€[« ,«i]. 

These results can be easily obtained by using the triangle inequality, the classical comparison between series 
and integrals, and the properties in (|6.2(l of the smoothing operators Sg . 



(7.3) 



Lemma 7.2. If 4 < sq < a and a > 5, then, for the modified state l/ ± >"+2 } we have 
\\V^ n+1 - - Sg n V ±M \\ s , T < C56 s n +1 ~ a for any s e [4, Sl + 2]. 



Lemma 7.3. Let sq > 5 cmd a > so + 3. For fc < n — 1, we /icwe 



where 



and 



L 2 (s) 



L 3 (s) 



\e±, k \\s,T < cs 2 e 



e±,k\U,T < C6 2 6 



for s € [8o,ai - 4], 
for s e [s , 8i - 3], 



Mh'&cit) < G5 2 e L k {s) A k for s e [s - Mi -4], 



max((s + 2 — a) + +So — ol — 1, s + So + 4 — 2a) 
max(so — a, 2(so — a) + 4) if a = s + 4, 
so + 2 — a if a = s + 2, 

max((s + 2 — a) + + 2(s — a), s + sq + 2 — 2a) 
so — a — 1 if a = s + 3, 
Sq — a if a = s + 2. 



if a^s + 2,s + 4, 



if a^ s + 2,s + 3, 



max((s + 3 — a) + + 2(sq — a), s ■ 
so — a — 1 if a = s + 4, 
sq — a if a = s + 3. 



so + 3 — 2a) if a/s|3,s + 4, 



(7.4) 



10 



Denote the accumulated errors by 

n— 1 n—1 n — 1 

fc=0 fc=0 fe=0 

Then, as a corollary of Lemma 7.3, we have 

Lemma 7.4. Let sq > 5, a > so + 3, and Si < 2a — sq — 1. TTiera 

' |S±,n|| s ,T < Co" 2 0„ /or s G [s , si -4], 

|^±,n||.,T < S 2 for se [s , si -3], (7.6) 
J|^n||H«(Mir) <CP9 n forse [*o-l,*i-4]. 

Lemma 7.5. For any so > 5, a > sq + 5, and si £ [a + 5, 2a — so — 1], we have 

'||^ = lkr<CA n e-- a - 1 (ll/, ± lla 1 r + 0. 

lkn|lH-+i(60x) < C^Anfl'— 1 , (7.7) 

[ll/tflkr^OPAnflr*- 1 
/or aZZ s > so . 

The proofs of Lemmas 7.2-7.3 and 7.5 will be given in Section 8. With these lemmas, we can now prove 
the following key result for the convergence of the iteration scheme. 

Proposition 7.1. For the solution sequence (5V ±k ,5^ ±,k ,S(j) k ) given by (|tj.21|) and (|6.34[l . we /lave 

'||(*V ± .*,«* ± .*)||. ) r + 11^*11^-1(6!^) <^- a_lA * /orsG [s ,si], 
|| J C(7±.*+i J $±.fc+i)7 ± ^ 1 -/±||. ) r<2^°- 1 /or sG [ So , si -4], (7.8) 
i ||B(7 + '*+ 1 > 7-'*+ 1 > ^+ 1 )|| ir .- 1 (MVr) <^r X /orse [s ,si-2] 
/or any k > 0, where = 9k+i — Ok- 

Proof: Estimate (|7.8|) is proved by induction on k > 0. 

Step 1. Verification of (|T%|) /or fc = 0. We first notice from (|57j) - [(5~%jl that (?7±,*±) satisfies the 
Rankine-Hugoniot conditions, and V r± '° = $ ± '° = imply that V ± '? = 0. Thus, SV^' satisfies the 
following problem: 

'L'^^SV^ = SeJt mSl T , 

S;^ a) V + ' O »^"' O ,*0°)=O onM?,, (7.!)) 
/^ ± °°|*<o = 0. 
Applying Theorem 4.2 to problem (|7.9|l . we have 

||W ±,0 |L,T < C^ s - a - 1)+ ||/ a ± || a+1)T for all s S [s ,si]. (7.10) 
Similarly, from i|(i.34fl . we obtain 

\\S^°\\ s ,t < C (\\ht\\ s ,T + \\U?\\.,t\\HU,t) , (7.11) 



where h* = SV 2 '° - d^^ol^ '° - d Xa ^SV^° satisfies 

ll^lkr < c\\sv±>% tT < c et°" x)+ \\ft\U + i, T . (7.i2) 

Thus, we have 

||^ ±,0 ||.,r<C 1 ^- a - 1)+ ||/ a ± || Q+1 , T . (7.13) 

Since 

5V ± = SV^ - S^ ?* 1 ^ 

we use (fTTTUf t and lf7~T3jl to find 



WSV^IUt < C 2 et a ~ 1)+ \\ft\Ui,T. (7.14) 
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From l|7rT3)) - lJTTljl . we deduce 

IKd^' ,^' )!!^ <^- a_1 Ao for any s€ [s , Sl ], (7.15) 

provided that \\f^ \\cc+i,t/S is small. 
Obviously, one has 

^±,^±,1)^,1 _ f ± h T < || (5flo _ 7)/ ±| kT + || e±)0 || SiTj (7.16) 

where 

e±,o - L(Ut + V^\9i + ^)(Ut + V^ 1 ) - L(Ut^t)Ut 

satisfies 

||e±,o|U,r<C^- a - 1)+ ||/±|| a+ i lT for se [s 0) si-2]. (7.17) 
Using the property of Sg , we have 

ll(Seo--0/a ± ll.,T < C6 s Q - a - l \\f±\\ a+hT for s<a + l, 

and 

\\(Se -I)f^\\s,T<Ce s - Sl \\ft\\s u T for.se [a + 1, 81 ]. 
Thus, from l|7.16|l . we obtain 

\\C(V ± >\$ ± ' 1 )V ± > 1 -f?\\ a ,T< 280' 1 - a - 1 forsG [s ,si-2], (7.18) 

when ||/^||a+i,r/5 is small and 9q > is properly large. 
Similarly, we can obtain 

\\B(V + '\V-'\ ( p 1 )\\ H ,- HmT) <^- - 1 for S G[ So ,si-2]. 

Sfep i?. Suppose that (I7.8|l ZioZcfe /or k < n — 1, we now; verify H7.8J1 /or fc = n. First, we note that, 
when a > so + 3, 

|| (t£ + V±.«+*,*± + S flB $ ± ' n )|| ao+3 , T 

< ||(t^=, *±)|U +a ,r + HV^'* 1 - 1 -* - ^^"lU+a,!. + IK^V^-", S , e„* ±,t *)]U+2,r (7-19) 

< C(5(l + 6>^ +3 - Q + e^ o+2 - a) +) < C5, 

by using assumption (|7.1(l and Lemmas 7.1-7.2. Applying Theorem 4.2 to problem H6.21J) . we have 
II^'IUt < C(||/±||,,t + \\9n\\ H * + i(bn T ) + l|coe/|| S;T (||/±|| SQ , T + \\9n\\ H °o +H m T ))) 

< C(ll/n IU,T + ||fl»||ff.+i(bn a .) + Wt + V^+K^ + S 9n $ ± '") IU+2,t(|| ft IL ,T + ||<7»|| ff-o + i (6n T ) )) ■ 

(7.20) 

On the other hand, similar to Ij7.19|l . from assumption l|7.1(l and Lemmas 7.1-7.2, we have 

\\(Uf + V^+i^t + 5 e „$ ± '™)|| s+2 , T < C5(l + 0°+ 3 - a + 0i s+2 - a) +), (7.21) 

which implies 

||({7± + V^+K^t + 5 e „$ ± '")|U +2 ,T(||/,tl| S0 ,T + ||5„||^o +1( 6n T )) < Ce^AniWftW^T + 5 2 ) (7.22) 

for all s £ [sg, si] by using a > sq + 3 and l|7.7ll . 
Thus, from [|7.20JI . we conclude 

||<5U ± ^||,,T<^r Q_1 A„(||/ Q ± |U : T + ( 5 2 ) forse [« ,*i]. (7.23) 
For problem (|6.34|) . we can easily obtain the following estimate: 

||5$ ± ' n |U,T < C(||fe±||,, T + 11^ + ^y^lU.rll^lUo.T) for any s > s , (7.24) 
where C > depends only on \\U^ + Se n V ± ' n \\ So ,T, and 

ht = ht + 8Vt' n - d X2 (*± + So***'")™?'" - d X3 (^t + Se n ^ n Wt n - (7-25) 
Obviously, we have 

H^lkT < \K\\ S ,T + (1 + ||*± + & n 3 ±,n |k,T)||^' n ||.,T + ||*± + S flB $ ± '"|| a+1 ,T||^' n ||.o,T, 
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which implies 

\\ht\\.,T < CA n ^-«- 1 (||/ a ± ||a,T + S 2 ) (7.26) 
by using Ij7.23|l and Lemmas 7.1 and 7.5. 
Substituting into lf7~2^|l . we find 

||<5$ ± ^|| s ,T<CA„ftr«- 1 (||/ a ± || aiT + ( 5 2 ) for alls G [s , si]. (7.27) 



Together (l7~2li|) with {?"2T|l . we obtain (EHJ) for ( < 51/ ± '™, <5$ ± '™) by using 



and choosing || q ,t/<5 and S > to be small. 

Finally, we verify the other inequalities in l|7.8|l . First, from (|tj.21fl - l|tj.22fl . we have 

£(^±,n $±,n )V ±,n _ f ± = _ j )f ± + (/ _ Sg^E^l + e± „_ 1 . ( 7 . 28 ) 

From Lemma 7.4, we have 

||(J - Sfl„_ l )£ t , Tl -i||.,T < C^ilSi.n-ills.T < CT,r 5+1 <5 2 < C^-"-^ 2 (7.29) 

by choosing s = a + 2. 

From Lemma 7.3, we obtain 

||e±,„_i|L,T < C9^ {s) 5 2 < C9 s - a ^5 2 for all s G [* ,«i - 4], (7.30) 

by using a > sq + 5. 

When s < a + 1 , we have 

\\(Se n ^ -I)fth,T < C0'- a - 1 \\ft\\a+i,T\ (7.31) 
while s G (a + 1, Si — 4], we have 

110%,.-! -J) At < Ilft n _ 1 / a ± ||.,r + 11/^11^ < ce^-'WftlU^T + s (7.32) 

by using l|6.2[) and l|7.1|l . 

Substitution fT^ - fTS^ into fT^ yields 

|| jC(v ±,« j$ ±, n)v ±,n _ /a± || s T < 2 ^«-«-i for s G [s 0) Sl - 4], 

provided that ||/^||si-4,t < S and ||/^ ||q+i,t/<5 is small. 
Similarly, we have 

\\B{V + > n ,V- n ,n Hs - HmT) <56 s - a -\ (7.33) 

Thus, we obtain (|T.8|> for k = n. 

Convergence of the Iteration Scheme: From Proposition 7.1, we have 

^(||(^ ± ' n ,<J$ ± - n )|| Q ,T + ||^ n || ff «-i (6nT) ) <oo, (7.34) 

which implies that there exists (V ± , $ ± ) € B a (Q T ) with G i2" a_1 (6ft T ) such that 

J(V ± '",$ ± ' n ) — > (V*,^) in S a (fi T ) x B a (Q T ), 
U n — ► </> in ff a - 1 (6fi T ). 

Thus, we conclude 



(7.35) 



Theorem 7.1. Let a > 14 and «i G [a + 5, 2a - 9] . Let if) e H 2si+3 (M 2 ) andU^ -U ± B 2[si+2 \R 3 + ) 
satisfy the compatibility conditions of problem l|2.17|l - l|2.20|l up to order Si + 2, and let conditions (|2.3|I - H2.4|) 
and 1|7.1[) 6e satisfied. Then there exists a solution (V^,^^) G _B Q (f2T) with <fi G H a ~ 1 (bQT) to problem 
(IB~T2l . 

Then Theorem 2.1 (Main Theorem) in Section 2 directly follows from Theorem 7.1. 
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8. Error Estimates: Proofs of Lemmas 7.2-7.3 and 7.5 



In this section, we study the error estimates for the iteration scheme (|6.21|) (I6.22[) and i|ti.34[l - H(i.35l) to 
provide the proofs for Lemmas 7.2-7.3 and 7.5 under the assumption that (|7.8|) holds for all < k < n — 1. 
We start with the proof of Lemma 7.2. 



Proof of Lemma 7.2. Denote by 



£f' n :=E(V ± ' n ,$ ± < n ) 



and 



±,n 



as the extension of (B{V ± ' n , (f> n ))^ in (l T . 

By the definition in l|(j.25fl for V^' n+ ^ , we have 



= S e jt' n + [dt,Se n ]^' n + (d X2 (*t + Se n ^ n )Sg n V^ n - Sg n (d X2 (H>± + $ ± '«)V 3 ± ' n )) 
+ (^ 2 (^„$ ± '")C/ a ± 3 - Se n (d X2 $±>"U± 3 )) + (d X3 (Se n ^ n )Ut A ~ Se n (d X3 ^ n U^)) , 



v ±,n+i_ ±, n 



+ (d X2 (Se n ^ n )Ut fi ~ Se n (d X2 ^ n Ut fi )) + (d X3 (Se n ^ n )U^ 7 - S Bn (d X3 ^ n U^ 7 )) . 
On the other hand, we have 



( c±,n 



±,n-l 



±,n 



+a :C3 (5$ ± ^- 1 )(^ a ± 4 + ^'"- 1 ), 



±,n-l 



sv 2 

±,n-l\ 



±,n-l 



+ ^ 2 (^ ± '"- 1 )(f/ a ± 3 + y 3 ± ' n - 1 ) 



±,n-l 



±,n-l 



+d X3 (6^-i ){u ± 7 + v ±,n-i^ 



±,n-l\ 



which implies 

by using A n = 0(9~ 1 ), the inductive assumption for l|7.8|l . and Lemma 7.1. Thus we deduce 



Se n (st n ^2' n )\U,T<cse s n 



for any s > Sq. 



(8.1) 
(8.2) 



(8.3) 



(8.4) 



(8.5) 



(8.6) 



The discussion for the commutators in (|8.3I) follows an argument from |2] . We now analyze the third term 



of Vr 



±,n+i 



Sg n V^ z ' n given by (|8.3|) in detail. 



When s e [a + 1, Si + 2], we have 



<c(||9 X2 (*± + 5^^)11^ ||s en v; 



and 



<cer a 11^(^+^)11^11^ 



±,m 



|*± + $ ± ' n ||a+i,r||V 3 ± ' n || i0 o < C<5 2 ^ +1 -« 



by using l|7.1|l . the induction assumption for (|7.8|) . and Lemma 7.1, which implies 

\\d X2 (*t + Se^nSe^ - 5 e „(S X2 (*± + $ ± ' n )F 3 ± '")|| s , T < CT 2 ^ 

20 



(8.7) 



When s € [so,a], we have 

\\d xa (?t +5 e „$ ± ^)5 e „l/ 3 ± '" - Se n (d^t + ^)Vt' n )\\., T 

< 11(1- S e J(d X2 (*± + ^ n )V^ n )\kT + R 2 (*± + ^' n )(Sg n - W^IUt 
+ \\d X2 ((Se n -I)<S> ± < n )Se n V^ n \\ s , T 

(8.8) 

+cei-«\\v±- n \\ a A^ + * ± ' n IUi,r + \\d X2 ((Sg n - m^ n )S9 n vf' n \\ s , T 

< CS 2 9f+ 1 - a . 
Together (JSTJ) with it follows that 

R 3 (tf± + S^'™)^^'" - Sg n (d X2 (^ + $ ± '")^ 3 ± ' ri )|| s , T < C5 2 e s - a+1 for all a G [««,,*! + 2]. (8.9) 

Other commutators appeared in (|8.3|l satisfy the estimates similar to the above. Thus, we complete the 
proof. 

To show the error estimates given in Lemma 7.3, we first show several lemmas dealing with different type 
errors. 

Lemma 8.1. Let a > So + 1 > 6. For the quadratic errors, we have 

||e£ 1} JU.r ^ C5 2 6 L k l{s) A k forse [s -2, Sl -2], 

llegfcll-.r ^ CT 2 ^- 2 " +S0_2 A fc /or s e [s - Mi - 1], (8-10) 
lle^lltf'(MJT) <C5 2 9 s - 2a+Sa - 1 A k forse [s -2, Sl -2] 

for all k < n — 1, where 

L 1 (s) = max((s + 2 — a) + + 2(s — a) — 3, s + so — f — 2a, s + 2s — 3 — 3a, s — a — 3). (8-11) 

Proof. The proof is divided into three steps. 
5"iep 1. From the definition of e± \, we first have 

= / (i-r)£V +v± .* +r5V± .».^^^ 



From (|7.1|) and Lemma 7.1, we find 

sup (HL^ + + r^^Hwri,.^) + ||*^ + t> ± k + r6^ k \\ w ^ {nT) ) < CS. 

0<r<l 

On the other hand, we have 

ll^ ± ,* ± )((^ ± ' 1 ^ ± ' 1 ),(^ ± ' 2 ^ ± ' 2 ))ll^ 

<c(||(c/ ± ,vi/ ± )|| s+2iT ||(y ± -\$ ± - 1 )|| wl , o||(y ± '^$±.2 ) ||^ 

+ ik^ 1 ,^' 1 )!!^!!^ 2 ,^ 2 )!!^ + ik^ 1 ,^' 1 )!!^.-!!^' 2 ,^ 2 )!!^) 

Therefore, we obtain 

llegjkr <C( ( 5 2 ^ (S "- 1 - Q£) ^|(5 + ||(V- ± - fc ,«I> ± ' fe )|U- +2 , r H- ||(5^ ± ' fc ,<5cI>±^-)|| s+2 ^) 
+ 5e s k °- 1 - a A k \\(6V ± > k ,6$ ± > k )\\ s+2tT ), 

which implies 

ll e ±*IU.r ^ C5 2 6 k {s) A k when s + 2 £ a, s< Sl -2, 
\\eQ k \\ a ,T < CS 2 e k n ^ is °- a - 3 ' 2so - 2 - 2a) A k when s + 2 = a, 
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(8-12) 
(8.13) 

(8-14) 
(8.15) 



by Lemma 7.1, the inductive assumption for 1)7. 8|) . and = O^." 1 ), where 

L(s) = max((,s + 2 — a) + + 2(sq — a) — 3, s + So — 1 — 2a, s + 2so — 3 — 3a). 

Thus, we conclude the first result in l|8.10[l . 

Step 2. Obviously, from the definition of e£ \ in l|t).29[l . we have 

\\e%\\ s , T <C (11^*11^11^*11^1,0. + \\8V^ k \\ L ^\\8^ k \\ s+l , T ) , (8.16) 

which implies the second result in (|8.10() . 

Step 3. Since ej_!\,Ui=o = (&k)t , we obtain that (efc)f satisfy the estimate given in i|8.10[l . Moreover, 
from the definition of , we have 



(ef ) ) 2 ± = -d^W^ - B^W?'", (gW) 3 = \{\5V^ - |^< fe | 2 

which implies 

fll(4 1) ) 2 ± ll^( b o r) < C(\\8^\\ Hs+ , {mT) \\8V^ k \\ L ^ + \\6<P\\ W i,~ (Mt) \\W±'H h . (UIt) ) 

lll(4 1) )3||^(60 T ) < C\\8V±*\\ L ~ 11^*11^(6^), 

yielding the last estimate of l|8.10|l by the inductive assumption for l|7.8|l . 
For the first substitution errors, we have 



where 



Proof. The proof is divided into three steps. 

(2) 

Step 1. From the definition of e^. k , we have 



-,(2) 



T 



L" , 



e ±.k - I ^ (!7 Q ± +S( ))c y±. fc +r(l-S ( ,Jy±. fe ;*J+S ( , fc *±- fc +T(l~S ( , fc )*±. fc ) 

x ((6V ± ' k ,5$ ± ' k ), ((1 - Se h )V^ k , (1 - S 9k )^ k ))dT. 



o 



As in l|8.13[) . from the inductive assumption for l|7.8l) . we find 



(8.18) 



Lemma 8.2. Let a > s + 1 > 6. Then, for all k < n — 1, we have 

'||ei 2) J s ,T<C<5 2 ^ 2(s) A fe / ors€ [ So _ 2 , Sl -2], 

tu\\s,T < C8 2 6 s - 2a+s °A k forse [s - l,s, - 1], (8.19) 



l4 2) ll^(6Qr) < C8 2 9l- 2a+s ° +1 A k forse [s - 2, Sl - 2], 



max(5 — a, sq — 1 — a) /or a = s + 2. 



(8.21) 



sup + Sg k V ± ' k + t(1 — Se k )V ± ' k \\w 1 ,°°(n T ) < C<5, (8.22) 

0<r<l 

and 



sup \\*$+S 0h &> k + T(l-S gh )$ ± > k \\ w i.~ (aT )<C6. (8.23) 

0<T<1 

Then we use ()8.14j) in 1)8.21)1 to obtain 

x (||({7±,*i>±)|U 2 , T + WSedV^^^h^T + ll(i - Se k )(y ± > k ) $±*)\\ a+2lT ) 

+ \\(SV^ k ,8^ ± ' k )\\ s+2 , T \\(l-Sg k )(V^ k ,^ ± ' k )\\ w i^ 

+ ||(^ ± ' fe ,^ ± ' fc )||K,i,»||(l--5 9 J(F ± » fe ,$ ± ' fc )|| s+2 ,r). 

Using the properties of the smoothing operators, the inductive assumption for (|7.8|l . and Lemma 7.1 in 
(|8.24(l . we conclude the first estimate in 1)8.19(1 when s < si — 2. 



Step 2. From the definition of e^\, in l lOty . we have 

\\e%h,T <chsV^' k \\ a A^^^ " ftJ^^IU- + ||(1 - Se^h+i^W^U- 

\ (8-25) 

+ ||(1 - S ek )V^ k \\ s , T \\V {x2 , X3) 5^ k \\ L ^ + \\S^% +1 , T \\(1 - SjJV^Hz-), 

which implies the second result in l|8.19|l when s < Si — 1 by the inductive assumption for (|7.8|) and Lemma 
7.1 in (H2S|)- 

Sfep 3. Noting that 

/-=(2)s± _ -(2) I 

\ e k )i ~ e ±,fcUi=o 

from the above discussion, we conclude that (&^)f satisfy the estimate given in (|8.19|l when s < s% — 2. 
From the definition of e k , we have 

'(4 2) ) 2 ± = d X2 ((So k lW)5V±- k + d X3 ((S 8k - I)<t> k )5V^ k + (So k - I)V^ k d X2 {5cp k ) 

+{Se k -I)V^ k d X3 {5cj> k ), (8.26) 
.(4 2) )3 = (I- Se k )V+- k 5V+- k -(I- Se k )V H > k 5V H ' k , 
which implies 

'iKsf^lk-CMV) < C(||*V ± . fc || H . (6nT) ||(l - S 9k )cf> k \\w^ + 11(1 - S 0h )<ffi\\ B . +1[mT) \\SV ± >''\\ L - 

+||(1 - So h )V±' k \\ H . iwlT) \\64fi\\ w i,~ + ||<ty fc || ff . +1(MlT) ||(l - SojV^Uao), 

ll(4 2) )3ll^(^ r) < (7(11^^11^(^)11(1 - Se k )V±' k \\ L ~ + ||(1 - Se k )V^ k \\ HHmT) \\SV^ k \\ L ^). 

(8.27) 

Using the inductive assumption for l|7.8|l and Lemma 7.1 in (|8.27|) . we conclude the last result in (|8.19|) 
when s < Si — 2. 

The representations of the second substitution errors {e±\, e^f' ) ) are similar to those of (e^L e k 2 ^). Thus, 
using Lemma 7.2 and the same argument as the proof of Lemma 8.2, we conclude 



Lemma 8.3. Let a > sq + 1 > 6. For the second substitution errors, we have 

( ^fjs,T<CS^ L k 



i 3) J s ,T < CS 2 6^ s) A k for se[s - 2, fll - 2], 
\\ef\\ H ^ b n T ) < CS 2 e s k +s « +2 - 2a A k for s e [s ~ 2, fll - 2] 



(8.28) 



for all k < n — 1, where L^^s) = max(s + sq + 2 — 2a, s + 2sq + 3 — 3a, (s + 2 — a) + + 2(so — a) + 1). 
Lemma 8.4. Let sq > 5 and a > so + 3. For the last errors, we have 

' \e ( ±\\\ s ,T < C5 2 9^ (s) A k for s e [* 0> *i - 4], 



f 1 



i 4) fe || s ,T < C5 2 9^ {s - l) A k for s e [*„, «i - 3], (8.29) 



Ul4 ll^»(6n r ) < C<5 2 ^ 5(s) A fc /or s e [«„ - l )Sl - 4] 
for all k < n — 1, where 

{max((s + 2 - a)+ + 2 + 2(s - a), s + 4 + s - 2a) if a/s|2,s + 4, 
max(so — a, 2(s — a) + 4) if a = s + 4, 
s + 2 - a) if a = s + 2, 



and 



max((s + 3 — a) + + 2(so — Qi), s + 3 + So — 2a) if a^s + 3,s + 4, 
^s(s) = s - a if a = s + 3, 

so — 1 — a) if a = s + 4. 

Proof. The proof is divided into three steps. 
Step 1. Set 

fl± =^ 1 (L(L/ a ± + U ± - fc+ =,vI/± + ^ fe <i> ± ^)([/ a ± +y ± ^)). (8.30) 
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Then we have 

\\Ri\U,T <\\L(Ut + V^+Kvt + Se k ^ k )(Ut + V^+i) 
- L{Ut + v±- k , *± + ^ k )(Ut + v^w^t 
+ \\£(V ± < k ,<S> ± > k )V ± > k -f±\\ s+2 , T , 
which implies that, for all s € [sq — 2, si — 4], 

||i£||.,T <C(||y ± ' fe +^ - ^ ± ' fc || wl ,c.(||L^ + K ± '* + »|| B+a ,T + + Se k ^%+4,T) 

+ \\tlt + V^Ws+^riWV^^ - V^Wl- + ||(7- Se h )^' k \\ w ^) 

+ \\Ut + V^ k \\ w ^(\\V^ k+ i - V ± '*||. + a,T + - Sg k )^ k \\ s+ ^ T ) 

'CT 2 ^ +4 ~ Q)++S0+1 ~ Q + CT^+ 5 - Q for s + 4 ^ a, 
C ( j2 6 ,so+2-a + C5 q for s + 4 = a. 



(8.31) 



< 

Thus we find that 



J 4 ) 



RtS^ ± ' k 



satisfy 

He^lU^ < C(||^|| S0 _ 2 , r (^ _1 ~ Q ^ fc + ^- 1 —A fe ( ( 5 + ^ S+2 — )+ )) + < 5^- 1 — ^H^IU^), (8.32) 
when s + 2 a, and 

He^lU^ < Cdl^lU^a^C^- 1 -"^,. + ^^-^"^feCJ + ^log^fc)) + ^^-^"^fcH^IU^), (8.33) 
when s + 2 = a, which yields the first estimate in H8.29J) for any s € [sq, si — 4], provided a > sq + 3 by using 

Step 2. Set 

R\ = B{ut + V^+Ki^a + SeJ k ). (8.34) 

Then we have 

\\R h k \\ H >(m T ) < \\B(Ut + V±< k+ i,i> a + S 9k 4> k ) ~ B{Ut + V^ k ^ a + <P k )\\m { m T ) + , fc )|| ff . (6nT ), 

which implies 

\\(Rk)i\\H*(bQ T ) 

< c(\\(Se„ - l)0 fe ||^+i(M M + \\vt' k+h - V^ k \\ H s (bnT) + \\(S 6k - 1)V ±m \\lA\^ + Se k <t> k \\ H ^ ( m T ) 
+ \\(S 6k - I)V^ k \\ H s {bnT) U a + SeJ k \\ w ^ + \\(S 6k - l)cf> k \\ Hs+HbnT) \\ut + V^ k \\ L ^ 
+ \\(S 9k - I^Ww^^WU^ + V^Wh,^ + HB^.^IU-cmv)). 

(8.35) 

Therefore, we have the estimate 

f C( J^»«((»+2-a) + +^-a,«+2-a) if a ^ s + 1, s + 2 , 

||(-Rfe)fll^(6fi r ) < | C86 k ifa = s + l, (8.36) 
[CS if a = s + 2 

for all s G [so — 1, si — 2]. Thus we find that 

[ k h " d Xl (^ + So k ^ k ) Xl =o^ 

satisfy 

IKef^fH^^) < C(||(^)f|| ff4(bf2T) (^r Q A fc + ( 5 2 ^" 1 - Q A fc ^ +3 - a)+ ) 
+^- 1 -«A fe ||(^)f||^ +2(fc0r) ), 
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(8.37) 



when s + 3 =/= a, and 
\\(ef)t\\ H s {mT) <C(\\(Rl)f^ 

(8.38) 

when s + 3 = a, which yields the estimate of (ej^)i given in (|8.29|l for any s 6 [s — l;Si — 4] by using 



Estimate (|8.29|) for the other components of el 4 ' can be proved by the same argument as above. 

Step 3. Noting that the restriction of e±\ on {x\ = 0} is the same as we conclude the estimate 

of e^ fc in (jQ9"|) for any s e [s , si - 3]. 

Proof of Lemma 7.3: Summarizing all the results in Lemmas 8.1-8.4, we obtain the estimates given 
in Lemma 7.3. 

Proof of Lemma 7.5. From the definitions of (f^,g n ,h^) given in 1)6. 22|) and l|tj.35fl . respectively, we 
have 



f ± = 

J n 


(Se n 


- ^e„-i)/a _ (Se„ - 


S0„-i)E±,n-l — S$ n e± t n-li 




9n = 




n — Se n -i)E n -i — Sg r 


c-n— 1 •> 


(8.39) 


h$ = 


(Se n 


-i — Sg n )E± :n ^i — Se 


„e±,n-l- 





Using the properties of the smoothing operators, we have 

'\\{Se n - S en _J/±|| s> T < Cd^AJfthr for 5 > 0, 
\(Se n - S^.JK^-iIUt < C^- s - 1 A n ||^±, n _i||g >T < CS 2 s n - s A n for s e [s , 8l - 4], (8.40) 
l^e^-xl^T < C^- s " )+ ||e ±>n _i|U-, T < C<5 2 ^ s " S)++Ll(S) A„ for s e [«„, *i - 4]. 

Estimate (|8.40fh can be represented as the following three cases. 
Case 1: s = sq. For this case, if we choose s S [sq 7 a — 5], then 

L 1 (s) = max(so — a — l,s + So+4 — 2a) = so — a — 1, 

which follows from (|8.4U|h that 

||S , fl B e±, n _i|| i0 ,T < C5 2 eir a ^^n- (8.41) 
Case 2: s = sq + 1 or s = s$ + 2. For these two cases, if we choose s = a — 4 > sq + 1, then 

L 1 (s) = max(so — a, 2(sq — a) + 4) = so — a, 

which implies 

||5fl B e ± , n _i||. lT < C6 2 6 s - a - 1 A n . (8.42) 
Case 3: s > so + 3. For this case, if we choose s = a — 2 > sq, then 

(s - s) + + L 1 (s) = (s + 2 - a)+ +s + 2- a<s-o;-l 
by using a > so + 5, which follows from (|8 .40^1 ^ that 

||Sfl B e±,„-i|L,T < CS 2 e^ a - 1 A n . (8.43) 

On the other hand, we have 

\\(S 6n - Sfl n _ 1 )^±,„_i||.,T < CS^r^An (8.44) 
by setting s = a + 1 in (|8.40() 9 if s% > a + 5. 
Together l|OT )l -l|Og )l with jO^I , it follows 

||/*|| s ,T<CTr Q - 1 A„( ( 5 2 + ||/±|| Q , T ) for all s> s . (8.45) 

Similarly, we have 

fll(5e„ - S e „_ 1 )E„_ 1 || K3+1(mT) < Ce s - s A n \\E n . x \\ H -s {mT) < C8 2 e s + 1 -~ s A n for S e [s - Mi - 4], 
\||ft n e n _i|| ir . + i (MlT) < C9i s+1 - S)+ \\e n ^\\ HHbnT) < C9i s+1 - S)++L3(S) 5 2 A n for S G [s -1,8!- 4]. 

(8.46) 

Estimate l|8.46|) o can be represented as the following three cases. 
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Case 1: s = sq. For this case, if we choose s — a — 4, then L 3 (a — 4) = sq — a — 1, which follows from 
(EUlb that \\SeJ n -i\\ff. 0+ i {b n T ) < CPO?-"- 1 ^ 

Case 2: sq + 1 < s < a — 3. For this case, if we choose s = a — 3, then L 3 (a — 3) = sq — a, which implies 
\\S 0n e a -i\\ H . + i (lbaT) < C5 2 s n °- a A n < C6 2 s n ~ a - 1 A n ; 

Case 3: s > a — 2. For this case, if we choose s = a — 5 > So, then £ 3 (s) = sq — 2 — a, which implies 
||S'e n e„_i||H-+ 1 (Mi T ) < C«5 2 0i^ s)++i3(5) A„ < C<5 2 ^ a ~ 1 A n . 

In summary, we obtain 

||Sfl B e n _i|| ff .+i (Mll o < C5 2 6 s - a - 1 A n for all s > s . 

From the assumption si > a + 6, it is possible to let s = a + 2 in I8.46|) i . which yields 

\\9n\\H°+i(bn T ) < CS 2 e s n - a - 1 A n for all a > a , 

by using and in fSSh- 

Furthermore, we have 

ll(5fl B - 5e B _ 1 ) J B±,n-i||. J T < C^- s - 1 A„p ±in _ 1 || s - r2 < C8H*- s - l A n for ,5 e [s 0)Sl 
||5 fln e ± , n _i||,,T < Ce4'~ a)+ ||e±, n _i||j,T < C0i s "" )++i2(s) <5 2 A n for s £ [s ,«i - 3]. 

Estimate l|8.49|l o can be represented as the following three cases. 

Case 1: s — sq. For this case, if we choose s = a — 3, then L 2 (s) = so — a — 1, which follows from l|8.49[) o 
that ||S , e„e±,„_i||. 0>T < C5 2 ^- Q - 1 A„; 

Case 2: s = Sq + 1. For this case, if we choose s = a — 2, then i 2 (s) = s — a, which implies 

||Sfl B e±, Tl _i||.o+i J T < C8 2 e s n °- a A n - 

Case 3: s > s + 2. For this case, if we choose s = a — 4 > Sq, then i 2 (s) = s — 2 — a, which implies 
||S , e B e±,n-i||.,T < C5 2 9 s n - a - 1 A n . 

In summary, we obtain 

||S fl „e±,„_i|| s ,r < C5 2 9 s n - a - 1 A n for all s > s . (8.50) 
Letting s = a in Q8.49|h . we find 

\\h±\\ s , T <C8 2 e s - a - 1 A n for all s> s , (8.51) 

by using |83U|) in JH31Jll- 



(8.47) 
(8.48) 



-3], 



(8.49) 
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